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Q Abstract. Spectral features absorbed against some radio quasars exhibit ~ 50mJy variations, with the lines 
varying both relative to the continuum and, when several lines are present, even relative to one another. We point 
£ — ■ out that such variability can be expected as a consequence of refractive scintillation caused by the interstellar 

medium of our Galaxy. Scintillation can cause independent variations between closely-spaced spectral lines, and 
can even alter the line profile. The background source need not be compact to exhibit spectral variability. The 
variability can be used to infer the parsec to sub-parsec scale structure of the intervening absorbing material. We 
discuss the importance of scintillation relative to other possible origins of spectral variability. The present theory 
■ is applied to account for the variations observed in the Hi-absorbed quasar PKS 1127—145. 

. Key words. Radio lines: galaxies - Galaxies: structure - Scattering - Techniques: high angular resolution 

O 1. Introduction 

i-^ ' There are now two well-established cases of variability in Hi-absorption spectra measured along lines of sight to high 
^"i, redshift quasars (Wolfe, Davis & Briggs 1982; Kanekar & Chengalur 2001). Changes in the amplitudes of spectral 
Q ' features relative to one another and to the continuum are observed on time scales of days to weeks, and there is 
no evidence for shifts in the absorption frequencies with time. Such rapid changes in absorption spectra imply the 



^ ' presence of structure in the intervening absorbing material on scales far below those attainable with conventional 
imaging techniques. 

The variations could be due to inhomogeneities in the absorbing source coupled with changes in the line of sight 
, to the background quasar or effects due to interstellar scintillation. Spectral variability due to changes in the line of 
5-h ' sight occurs if the background source moves at high angular speed with respect to the absorbing system. The time 
scale of spectral variations is determined by the angular speed of the background source relative to the angular scale 
of fluctuations in the intervening absorbing system. 

The interpretation of the variability in terms of interstellar scintillation is more troublesome. No changes in the 
absorption strength relative to the continuum are expected from the simplest scintillation model, in which radiation 
from a point-like source is absorbed by foreground material. The next simplest model consists of a background source 
composed of several components compact enough to scintillate, with different spectral features associated with each 
source component. Variations in spectral features relative to the continuum are expected in this model, but appear 
unable to quantitatively account for the actual variations observed in these sources. For PKS 1127—145, Kanekar 
& Chengalur (2001) rule out a source model consisting of one point-like variable component and one component of 
constant flux density; the model predicts that when the source exhibits the same continuum flux density it would also 
exhibit the same absorption spectrum, which is not observed. Moreover, a source model composed of unrealistically 
many compact background components is required to account for the large number (> 4) of independently varying 
spectral features that are actually observed. An alternate scintillation model proposed to explain HI variability in 
pulsars (Gwinn 2001) can also be ruled out, as it requires the source to exhibit (diffractive) scintillation in the 
absorbing medium, which is not observed here. 

However, existing models do not take into account effects related to the finite angular diameter of the background 
source. The variability characteristics of the scintillation depend in detail on the source image structure. This is 
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intervening absorption system, z abs Scattering medium in our Galaxy 

Fig. 1. The geometry of the absorbing system. 



important because, even for the most compact quasars, the presence of any sub-structure in the absorbing medium 
influences the apparent image of the absorbed source. Thus the source image that scintillation responds to differs 
between the continuum and any spectral feature. The implication of this simple point has not been appreciated in 
the discussion of spectral variability, and it has a profound influence on the variability of spectral features due to 
interstellar scintillation. 

In this paper we compute the effect of scintillation on absorption spectra, taking into account the effect of inho- 
mogeneous structure imposed by the absorbing material across an extended background source. The theory detailing 
the response of scintillation to image structure is presented in Sect. [21 In Sect. [21 the theory is applied to compute 
the variability properties of some commonly encountered distributions of absorbing material. The simplest model of 
absorption against a compact background source is used to account for the variations observed in PKS 1127— 145 in 
Sect. 01 The implications of scintillation for the interpretation of absorption spectra are discussed in Sect. 03 and the 
conclusions are presented in Sect. El Throughout the text, we focus on systems observed in absorption. However, the 
theory applies to equally to spectral lines observed in emission without a background source and we indicate how the 
results should be interpreted for such systems. 

2. Theory 

Consider the radiation from a background radio source whose spectrum is altered by absorption lines due to a fore- 
ground system, as shown in Fig. 1. As it enters our Galaxy, scattering by the diffuse ionized component of the interstellar 
medium (ISM) induces variability, or interstellar scintillation, in the intensity of the radiation. The properties of the 
variability depend sensitively on the angular brightness distribution toward the background source, which is influenced 
by the angular distribution of the absorbing material (see Fig. 2). In this section we describe how a given distribution 
of absorbing material determines the characteristics of the spectral variability due to scintillation. 

We restrict the discussion to the effects of scattering in our Galaxy. Scintillation in the ionized medium of the 
absorption system may, in principle, also induce spectral variability. However, because the distance between the 
scattering and absorbing medium in this case is small, the apparent size of the HI absorbing structure is too large for 
it to be subject to the effects of scintillation. 

2.1. The scattering model 

For present purposes it is convenient to idealise the interstellar scattering material as being confined to a thin phase 
screen located at a distance L from the observer. This assumption is commonly used in astrophysical situations 
(Prokhorov et al. 1975; Rumsey 1975, Mercier 1962, Salpeter 1967), and a treatment incorporating the extended 
nature of the scattering medium introduces no additional effects relevant to the present discussion. 

It is necessary to characterise the phase fluctuations responsible for the scattering. Observations along many lines 
of the sight through the Galactic interstellar medium indicate that the power spectrum of density fluctuations in the 
ISM often follows a power law with index —j3 close to the Kolmogorov value of f3 — 11/3 (e.g. Armstrong, Rickett 
& Spangler 1995). The amplitude of the spectrum depends on the particular line of sight through the Galaxy, and 
is given by the scattering measure SM = J dzCjf(z), where Cj^(z) denotes the level of turbulence as a function 
of distance along the ray path (see, e.g. Rickett 1977). The information on phase fluctuations in the thin screen 
approximation is succinctly embodied in the phase structure function, D^r) — ([cj)(r' — r) — 0(r')] 2 ), which takes the 
form D^(v) = (r /r^m) I for (3 < 4 provided that r exceeds the inner scale of the power law (Coles et al. 1987). The 
diffractive length, r^iff, is the scale at which the root mean square phase difference across the phase screen is unity, 
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and it is defined formally in terms of the scattering measure as 
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where the scattering measure for a medium of thickness AL is 
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Armed with a description of the phase fluctuations on the phase screen, one can describe the power spectrum of 
intensity fluctuations across the observer's plane for a point source. The source is assumed to be located at infinity 
- an excellent approximation for the extragalactic sources under consideration here - so that the wavefronts incident 
on the scattering medium may be considered planar. For scintillation in the regime of weak scattering and refractive 
scintillation in the regime of strong scattering the power spectrum of intensity fluctuations is (Codona & Frehlich 
1987; Coles et al. 1987) 
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Strong scattering occurs when rdiff < y/L/k (e.g. Narayan 1992). In this regime the spectrum in Eq. J3J is cut off by 
the exponential term at wavenumbers k = r^ss (L/k) -1 = rC; , where r re f is the scale of refractive scintillations. In 
the strong scattering regime the exponential term becomes important before the sine term oscillates rapidly with n. 

Observations of pulsars and intra-day variable radio sources located off the Galactic plane indicate that objects 
viewed through the ISM at frequencies below ~ 4 GHz are in the regime of strong scattering (Taylor & Cordes 1993; 
Walker 1998; Kedziora-Chudczer et al. 1997, Dennett-Thorpe & de Bruyn 1999). Thus observations of HI absorption 
lines, for instance, are always in the regime of strong scattering. The refractive scintillation that occurs in strong 
scattering is a broadband phenomenon, which is to say that a point source exhibits similar intensity variations over a 
bandwidth comparable to the observing frequency (i.e. over a range Av ~ v). 

Ultra-compact sources can also exhibit another manifestation of strong scattering, called diffractive scintillation. 
The contribution of diffractive scintillation is omitted from the present discussion. This narrowband phenomenon 
contributes at high wavenumbers k ~ r^ is » r~ f . This is neglected because terms that contribute at such high spatial 
frequencies are suppressed unless the source angular diameter is below r<afi / L which requires source sizes 9 < 1 //as at 
v ~ 1 GHz. Not even the most compact extragalactic sources, including intra-day variable sources, appear to satisfy 
this requirement (e.g. Dennison & Condon 1981; Kedziora-Chudczer, Macquart & Jauncey 2001). 



2.2. Variability of an extended source 

In practice no extragalactic radio source can be regarded as point-like, and the power spectrum given by Eq. J3J needs 
to be modified to take into account both the angular structure of the extragalactic source and its foreground absorbing 
structure. The power spectrum of intensity fluctuations across the observer's plane is the product of the point source 
scintillation spectrum with the power spectrum of the source angular brightness distribution (Little & Hewish 1966; 
Salpeter 1967; Rumsey 1975): 



Pi{k;v) 
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The quantity 14 pp (r; v) is the visibility of the source image at frequency v, or equivalently, wavenumber k = 2irv/c: 



y app (r; v) = J d 2 9I app (9; v) exp (-ikG ■ r) , (6) 

The apparent image is determined by both the angular brightness distribution of the background source, / src (0; v), and 
the angular distribution of the opacity, t{9\ v\ due to the intervening absorbing medium: I app (8; v) — / S rc(0; v) e~ T ( ' v \ 
The structure of the absorption material has a direct influence on the variability of its associated spectral line. The 
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(a) (b) (c) 

Fig. 2. The variability due to the continuum image of a source (a) is different relative to the image of the source at two 
spectral lines at frequencies, say, (b) f a bsi an d (c) ^ a bs2- Since the observed variability properties depend sensitively 
on fine-scale source structure, the intensity scintillations caused by scattering in the ISM of our Galaxy are different 
in all three cases. 

scintillation properties of a source can change sharply with frequency when the angular distribution of the optical depth 
t(9; v) varies with frequency. This is important if the absorption at each spectral line is associated with a distinct 
angular distribution of absorbing structure. The background source itself need not be compact for the apparent image 
to contain compact features, since the absorption can also introduce compact features into the apparent image. Thus 
the background source itself need not be compact for spectral features to exhibit scintillation, as shown in Sect. I3~T"T1 
and Sect . 13 . 2"T1 below . 

In practice an observer does not measure intensity fluctuations across a plane, but rather the temporal intensity 
fluctuations at a single point on this plane. The two are related under the assumption that the phase fluctuations are 
frozen into the scattering screen (e.g. Rumsey 1975). The spatial power spectrum of intensity fluctuations, Pi(k), is 
related to the auto-correlation of the intensity fluctuations by a Fourier transform 

Ci(t;u) = (A/(t»A/(t' + t; v)) = j d 2 ^* 1 "?^;^, (7) 

where v is the velocity of the interstellar medium transverse to the line of sight to the source. A typical value of 
v = 30kms _1 is derived from direct measurements using intra-day variable radio sources (Dennett-Thorpe & de 
Bruyn 2003; Bignall et al. 2003). 

The foregoing theory applies to the autocorrelation of the intensity fluctuations at a single frequency. However, as 
we are interested in comparing the intensities of spectral features relative to one another and to the continuum, we 
also compute the covariance between intensity fluctuations at two distinct frequencies. 

The frequency dependence of the covariance is potentially influenced by both its dependence on the scintillation 
power spectrum, _P sc j llt (K), and on the variation in image structure with frequency. However, the scintillations them- 
selves are broadband (e.g. Goodman & Narayan 1989), which is to say that a points source exhibits similar intensity 
fluctuations over a large range of observing frequencies Av ~ v. Variations in P S cint('«) with frequency are therefore 
negligible for present purposes, and it suffices to consider only the influence of variation in image structure with 
frequency. The cross power between observations at frequencies v and v' is (Little & Hewish 1966; Rickett 2000) 



Pj(k; v- v>) = P scint (K)V app f — ; K* pp ^ — ; u' j (8) 

This is the generalisation of Eq. JSJ , with the power spectrum of the image at frequency v replaced by the cross power 
spectrum of the images at v and v 1 . The corresponding temporal cross-correlation of the intensity fluctuations is 



d(t; v, v') = (AI(t'; u)AI(t' + t; u')) = J d 2 K e lK v *Pj(K; u, v'). (9) 
The auto- and cross-correlations given in Eqs. and JHJ relate the scintillation properties to the image structure. 



2.3. Assessment of variability relative to the continuum 

Spectral line variations can occur in sources whose continuum level also varies. Since only variations relative to the 
continuum indicate the presence of substructure in the absorbing material it is necessary to construct a measure of 
the variability that removes the contribution from continuum variability. 
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It is possible to separate spectral line variations from continuum fluctuations by computing how the flux density of 
an absorption line would vary due solely to changes in the continuum level. The fluctuations of a spectral line whose 
flux density varies proportional to the continuum level, A/ c jep(' / abs) are 

A/ dop (^ abs ) = e-AI(v) = ( SaIC ~ Saha ) AI(y), (10) 

\ >Jsrc / 

where r is the optical depth, AI{v) are the continuum flux density variations, S src is the mean source continuum flux 
density, and S^bs is the mean absorbed flux density at the frequency ^ a bs of the absorption line. 

Variations induced by scintillation that are independent of the continuum level might be mistaken for variations 
in the optical depth, r. The amplitude of spectral line fluctuations independent of the continuum fluctuations is given 
by ATjndept^abs) = A/(f a b s ) — A/<iep (^abs) • The auto-correlation of independent intensity fluctuations is then 

Cindcp(^abs) = {AI mc \ cv {t' , ^abs) AL m i ep (t + t, f a bs)) 

= Ci(t; z/abs, ^abs) - ( ^ S1C „ *^ abs ) [Ci(t; v ahs , v) + Ci(t; v, v ahs )] + [ ^ Slc ^ abs j Ci(t;v,v). (11) 

\ ^sre / \ '-'src / 

This expression subtracts out the continuum variability contribution to the variations observed in the spectral line. 
It is zero for any source whose spectral features vary identically to the continuum fluctuations, as would be the case 
if the image of the source at J/ a b s were identical to its continuum image, if the absorption structure were distributed 
homogeneously. 

A similar quantity assesses the covariance of independent variations between two spectral lines. For lines at fre- 
quencies ^absi and ^ a bs2 with absorbed flux densities Sabst and S a hs2 the cross-correlation is 

Cindep (^? ^absl 7 ^abs2) = (A/i nc }ep(^ , J y absl)A/i n dcp2(i' + t, f a bs)) 

sy /. \ , I S S rc Sabsl I ( '-'src >S'abs2 \ f, /. s 

= Cj(t, I/absl, ^absa) + o I o J C I\t,V,V) 

\ »^src / \ '-'src / 

Ssl '% Sahs2 ) d{t, u &hsl ,v) - ( S «°- S ** \ Cl (t, „, v ahs2 ). (12) 

'-'src / \ ^>src / 

This quantity removes fluctuations in the source continuum flux density when the covariance in the intensity fluctu- 
ations between two spectral features is computed. Again, it is identically zero if a spectral line and the continuum 
exhibit identical fractional flux density variations. 



2.3.1. Relevant angular scales 

Before applying the foregoing theory to specific absorption distributions, it is appropriate to consider the range of scales 
probed by scintillation. Scintillation is most sensitive to image structure on angular scales comparable to rc f = r rc f / L 
(viz. Eq. . The refractive angular scale is expressed in terms of standard scintillation parameters as follows 

/ r 2 \ 3/5 / j \ 3/5 / < xll/5 

^ = 2.15( TF ^_) (_) (_) mas , (13) 

where we have equated the thickness of the scattering material with the distance to the scattering screen, AL = L. 
For instance, scintillation of a spectral line due to neutral hydrogen at z = 0.1 (at A 0.232 m) probes structure on a 
spatial scale of ~ 2pc, for typical scattering parameters Cff = 10~ 4 m~ 20 / 3 and L = 0.5 kpc for a line of sight off the 
Galactic plane. The response of scintillation to structure on scales larger than # ro f decreases sharply. The modulation 
index of a source of size 8 scales as (6*/6' ro f)' 3 / 2_3 once its size exceeds 6* rc f . 



3. Absorption Models 

We now apply the foregoing results to two hypothetical distributions of absorbing material: one a population of clouds 
with fixed angular size, the other a medium whose power spectrum of absorption brightness fluctuations follows a 
power law. Absorption against both extended and compact sources is considered. 

The models are also applicable to variability from systems in emission. A distribution of material in emission would 
exhibit the same variability properties as the same distribution seen in absorption against an extended background 
source. This is because the power spectra of the angular brightness distributions in the two cases are identical, except 
for the contribution at zero spatial frequency (k = 0) due to the background source. As mentioned above, this 
contribution does not influence the variability. 
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3.1. Structure with a Characteristic Angular Scale 

3.1.1. Emission, or absorption against an extended background source 

Consider an extended background source whose spectrum is altered by a single intervening system. The intensity of 
the background source, Iq, is assumed constant over an angular scale much greater than the refractive scale $rcf, so 
that it does not scintillate. 

We consider a model in which the intervening system is comprised of two sets of clouds, with Ni clouds of size Q\ 
absorbing at centre frequency v\ and N 2 of size 9 2 absorbing at centre frequency v 2 . This is the simplest model which 
demonstrates that different spectral lines can vary relative to one another. The generalisation to an arbitrary number 
of spectral lines and cloud types is obvious, but introduces no new physics. 

To aid in the computation of scintillation effects the exponential of the optical depth is specified in the following 
idealised form 



exp[-r(0; v)) 



1 — A\ip{y — v\) exp 



1 - A 2 tp(v - v 2 ) E exp 



^ ' '9- Ad- 



2,3 J 



(14) 



where the angular offsets A9\i and A9 2 i describe the positions of the cloud centres on the sky. The absorption profile 
is described by the function tp(y), normalised such that the absorbed flux of a source of unit flux density is unity, 
J dv ip(v) — 1. We make the simplifying assumption that the line profiles at frequencies v\ and v 2 do not overlap. 
This ensures that the scintillations observed at frequency v\ are entirely due to one set of absorbing clouds, while the 
scintillations observed at v 2 are due to another set. The absorption strengths A\^ 2 determine the amplitude of the line, 
with < Ai t 2 < 1. The image visibility is then 

y app (r; v) = — hr^ - J o E A ^ n ex P — f- <p( v - ^) E ex P H** ■ Ae ^ > ( 15 ) 

where the two cloud types are indexed by 77. The first term represents the contribution from the background source, 
assumed to have a spatially constant flux density. This term only contributes at r = 0, and is henceforth neglected 
since it has no influence on the power spectrum of intensity fluctuations (see Eq. |@J). 

When the number of absorbing clouds is large one is justified in averaging over their positions. The power spectrum 
of the apparent image contains an average over the cloud positions (exp [— ikr ■ (A9 v ,i — A9 v j)]) with r] = [1, 2] and 
i.j = [1, . . . , N v }. The evaluation of this quantity has been treated in a similar context elsewhere (Macquart 2004), to 
which we refer the reader for more detail. 

It is convenient to separate this average into the N v self-terms and N„(N„ — 1) cross-terms. Each self-term yields a 
contribution of unity. The average over the cross terms proceeds by considering the clouds to be distributed randomly 
over an angular area fi, with p 2 (0i,9j)/£l 2 being the joint probability of finding the ith cloud at position 6i and 
the jth cloud at position 6j . We assume that the statistics of the cloud positions vary only weakly over the area in 
which the clouds are distributed, so that we may write the probability distribution as a function of object separation, 
AO = 6i — 9j only, and neglect its dependence on the average cloud location, © = (6{ + 9j)/2. Thus one has 

(exp [-ike ■ (Ad^ - A0„,i)]) = N v + Nr> ^ 2 ~ - f d 2 @d 2 A8p 2 {A0) exp[ikr ■ AO] 

NJNr, - 1) p 2 (ke) , s 

= n n ' (16) 

where p 2 is the Fourier transform of p 2 . 

The average power spectrum of the brightness distribution is 
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where S a i, ari — A v tt9 2 N v Iq are the flux densities absorbed by each population of clouds. 

In the simple case in which the clouds are sparsely distributed (i.e. N n /Q < 1), it is convenient to neglect the 
effect of the term involving p 2 . In fact, this term may always be ignored when the region over which the clouds are 
distributed greatly exceeds the refractive angular scale. In this case p 2 (r) falls to zero at sufficiently small r that it 
yields no contribution to the scintillation variations (cf. Eq. J7J). This is because the term involving^ only contains 
power on angular scales much larger than that probed by scintillation. 
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We now evaluate the auto-correlation of the temporal intensity fluctuations due to scintillation. We make the 
simplifying assumption that the angular scale of the absorbers is greater than or comparable to the refractive scale 
0i,2 > ^ref) which enables us to write the auto-correlation function in closed form (Coles et al. 1987). In this case the 
source power spectrum cuts off the integral in Eq. Q over k before the exponential term involving (see Eq. flf) 
becomes important and while the argument of the sine function is still small. Evaluation of the integral yields 

2 

-6 



(18) 
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v 2 t 2 
2' ' 2£ 2 02 



The variance of the intensity fluctuations is obtained by setting t = 0. 

In this case a numerical estimate for the root mean square flux density variation when j3 = 11/3 is 

1/2 

(19) 
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q 2 



— 1/2 

For a fixed absorbing flux density, the modulation index decreases as N n . This is because the observed intensity 
fluctuations may be regarded, in a statistical sense, as the scintillation pattern from independent absorbing clouds. 
The assumption that the absorbers are distributed homogeneously in a statistical sense over the area subtended 
by the background source ensures that the scintillation patterns due to individual absorbers are indepedent when 
averaged over the absorbing population. For a fixed number of absorbing clouds, Eq. (|19fl implies that large absorbing 
clouds contribute the most to intensity variations. The root mean square intensity fluctuation is proportional to 
N^ 1 ^ 2 9~"' 'Z 6 S'abs oc N^ 2 Q 5 / e . Thus the number of clouds of size 9 must rise faster than N„ oc (9~ 5 / 3 for decreasing 9 
for the contributions of small clouds to be equally important as those of the largest clouds. Otherwise, only the largest 
clouds in the absorbing medium contribute substantially to the observed spectral variability. 

We also derive the time scale of fluctuations, defined as the time scale at which the correlation function reaches 
1/e of its maximum value, from Eq. I|18|) 

L n / 



Wt„ = 36.9 —2- — — — - days. (20) 

' \1 mas/ V0.5kpc/ V30kms 1 J 

This time scale depends only on the size of the absorbers assuming 9 n > 9 re f. The ratio of fluctuation time scales of 
any two spectral lines is equal to the relative sizes of the absorbing objects, t Vl /t U2 — #i/02- 

The cross-correlation of the intensity fluctuations C; n dep(i, ^2) is determined by the cross power of the visibilities 
measured at v\ and i>2- This depends on the distribution of the positions of clouds of the two different sizes, 9\ and #2- 
A full description of the intensity cross-correlation involves the joint probability distribution of cloud centres between 
the two cloud types. However, a practical simplifying assumption is that the two types of clouds are distributed 
independently. If we assume that the cloud positions follow identical distribution functions, with the probability of 
finding a cloud i at position 6i being pi(6i)/Q, then the cross-correlation of the intensity fluctuations is proportional 
to 

^exp [zfcr • (0 M - (hj)]^ = ^|pi(fcr)| 2 . (21) 

Now if the clouds are distributed out to some maximum scale 6* m a X in the absorbing region, the quantity |pi(fcr)| 2 
declines rapidly for r > $max- I n anv realistic medium the absorbing clouds are distributed over an angular extent 
much larger than the angular scale probed by refractive scintillation (i.e. 9 max ^> rc f). Thus the visibility cross-power, 
Va.pp(KL/k, v\)V* pv {nL/k, v%), is cut off at wavenumbers k > {L9 raax )~ 1 . This wavenumber is much smaller than those 
probed by interstellar scintillation, n ~ r~ f (viz. Eq. (4)). This implies that the covariance of the intensity fluctuations, 
as described in Eqs. and (|12fl . is zero. 

In physical terms, the result in Eq. I|21(l shows that the cross-power of the visibilities at the two spectral lines 
exhibits structure only on angular scales comparable to the entire region over which the clouds are distributed. As this 
scale is much larger than that probed by interstellar scintillation, it does not respond to the image structure. This is 
the reason that the cross-correlation between spectral lines emanating from the two independently distributed sets of 
absorbing clouds is zero. 

In summary, we have shown that spectral lines due to absorbing clouds randomly distributed in front of a large 
background source can exhibit substantial variations due to interstellar scintillation. However, when multiple spectral 
lines are present the variations between the lines are not correlated under most conditions. 
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3.1.2. Absorption against a compact background source 

If the background source is itself compact enough to scintillate this has two important implications for the spectral 
variability: (i) the continuum level also varies, so that the assessment of spectral variability needs to be made relative 
to the varying continuum, and (ii) the radiation probes a only narrow cone of HI structure along the line of sight to the 
observer, so that only the absorption of relatively few clouds need be considered. To illustrate the effect of background 
source size on spectral variability we consider a background source with a Gaussian intensity distribution 



4rc(0) = h exp 



9 2 



(22) 



whose angular size is sufficiently small that it exhibits scintillation. The source background flux density is S src = IoTr9 2 rc . 
We consider absorption from the same intervening system of absorbing clouds as above, so the power spectrum and 
cross powers of the absorbed image are 



V app (r;v)V: pp (r;v) 
V app (r;v)V: pp (r;v') 
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(25) 



Variability from a collection of absorbers is treated in Appendix^] but in many practical cases it suffices to consider 
the intensity modulation of a background source with only a single absorber directly in front of it (A6 = 0). If both 
the absorber and background source are large, 9 srci 9\ ^ 9 rc { , we can use the same approximation made to derive Eq. 
I|18|) . The intensity auto- and cross-correlations are 



Cj(r;v) = 47r 2 r 2 A 2 SM (£\ 2^T U-^\ S 2 . c ^{L9 Brc f-\F 1 U 



-2A lip ^-^)(L9 CIOSS1 f- 6 
+A 2 1 <p(v-v 1 )(L9 cS f- 6 



■)2 
src 



A 3 



9 2 . 








2L 2 9 cS 2 



and 



G/(r;z/, v') = 4 7 r 2 r 2 A 2 SM(L/£ : ) 2 2 2 ^/ 2 r (3 - S 2 rc <p{v' 



-A (L9 C ioss) 



0-6 



9 2 

"src 



vi) (L9 SIC ) 



P-6 T? 

1-T1 



A 3 



p 1 r2 

2' ' 2L 2 9 2 



r 2 
2' ' 2L 2 9 2 



(26a) 



(26b) 



respectively where, 

src( 2 #l + ^src) 

) ■ 



9 2 

cross 



y eff — 7a 



9 2 

src 



■2 
src 



(27a) 
(27b) 



are effective angular diameters. The quantity 9, 



2 ross tends to 9 2 rc /2 when the background source angular diameter 
exceeds that of the absorber in the limit 9 STC ^> 6„, and tends to 9 2 IC in the opposite limit. Equation (|26a|l provides a 
means to compute the modulation index of the continuum fluctuations, line fluctuations and the "independent" line 
fluctuations. 

The time scale of the intensity cross correlation, set by # S rc and Across, is insensitive to the angular scales of the 
absorbers if they are smaller than the background source, since Across ~ 9 S rc/V^ is insensitive to 9\ in this regime. 
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3.1.3. Clumpy cloud distribution 

Equations 1(23(1 and 1)24(1 simplify when the positions of the clouds AO^j do not greatly exceed the background source 
size. In this case, the terms A<? 2 j in the exponentials of these two equations are negligible, and an average over random 
clouds centres yields 



(V app (r;v)V*(r;v)) = S, 



2 

src 



exp 



2„2/)2 1 



k 2 r 



77=1,2 



e 2 

src 



e 2 



N v exp 



k 2 r 2 9 2 rc 2 



and 

(V app (r;v)V: pp (r;v')) = S 2 IC exp 
The intensity auto- and cross-correlations are then respectively 

/*-« , « v 2 t 2 
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src 
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iFi\ 3-^,1, 



17=1,2 



7 c fF?7 



0-6 



e 2 

"src 



i^i 3 



ft 1 v 2 t 2 
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and 



where the quantity 



^ ef = 2^(/3-2)(2-?)rm(^ 



^diff 

'"F 



8-2/3 



(28) 



(29) 



(30a) 



(30b) 



(30c) 



approximates the variance that a point-like source of unit intensity would exhibit. The amplitude of the intensity 
variations exhibited by a spectral line, given by setting t = in Eq. (|30a|) . contains a component that is due to the 
size of the background source, and another component which scales linearly with the number of absorbing objects, 
N v . Note, however, that the amplitude of the cross-correlation between spectral frequencies v and v 1 is independent 
of the characteristics of the absorbers, and is merely equal to that which a continuum source would exhibit. 



3.2. A statistical description of emission or absorption fluctuations 

Instead of regarding the absorbing (or emitting) medium to be composed of discrete clouds, as discussed above, it is 
often more useful to consider the optical depth as a stochastic variable which is allowed to vary over the angular extent 
of the background source. In this section we derive the properties of spectral line variability caused by an absorbing 
medium whose spatial variations are specified only in a statistical sense. Spectral line variability in this model is caused 
only by spatial optical depth variations, which in turn may induce structure in the image of the absorbed source on 
sufficiently small scales to be subject to interstellar scintillation. The absorbing medium is not regarded as moving 
across the extent of the source itself, so no spectral variability is caused by relative motion of the absorbing material 
across the line of sight to the source. We assume that the optical depth variations associated with a given spectral line 
are confined to a thin screen a distance -D ac . s from the scattering medium in which the scintillations take place. 

In general terms one can describe the effect of the local medium on the incident electric field of the background 
source, w src (r), as follows 

u(v)=u SIC (r)e i ^ r \ (31) 

where the phase factor <f> is regarded as a complex number. The quantity Re(</>) is the phase delay caused by refraction, 
and Im((/>) = r/2 represents the effect of absorption (or amplification) by the medium 1 . For present purposes we are 
only interested in optical depth variations, so we henceforth set Re(0) = and consider only fluctuations in the 
quantity r\ = lm(0). 

Two specific assumptions are made concerning the statistical nature of the optical depth fluctuations. Firstly, it 
is assumed that the fluctuations are wide sense stationary, so that the statistical properties of r(r) do not depend on 

1 Note that Eq. l.'-iU reduces to the familiar result 7(r) = u(r)u*(r) = J src (r) e _T< - r '. 
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the position. Specifically, this assumption implies that quantities involving changes in optical depth between positions 
say, r' + r and r', depend only on their separation, r. 

The second assumption is that the statistical fluctuations in the optical depth follow a f distribution. We make 
this assumption in lieu of more detailed observational measurements of the distribution of optical depth fluctuations 
for any given spectral line. The assumption of Gaussian statistics has proven highly successful in treating fluctuations 
in the electron density variations in the ionized component in the interstellar medium of our Galaxy, which give rise 
to variations in Re(</>) (see, for example, Armstrong, Rickett & Spangler 1995). 

However, the model contains an implicit shortcoming, because normally distributed random variates are not con- 
strained to either positive or negative values. As such, it makes no distinction between the statistical properties of the 
spatial distribution of r depending on whether it is positive or negative. Nonetheless, this shortcoming is not serious 
if the optical depth variations are small compared to their mean value, so that most of the medium is seen chiefly in 
either absorption or emission. 

Recall from the discussion in Sect. 2.2 that the scintillation characteristics of the absorption (or emission) lines 
depend on the mean square visibility of the absorbed (or emitting) image. We derive a general result for this quantity 
in Appendix El and merely quote the result here: 



(|Kpp(r)| 2 ) = \V SYC (v)\ 2 exp 



-2(r> + -C T (0) - D T (r) 



(32) 



where V^ rc (r) is the visibility of the background source, C T (r) = ([r(r + r') — (r)][r(r) — (r)]) is the autocovariance of 
the optical depth fluctuations in the absorbing medium and D T (r) = 2[C T (0) — C T (r)] is its corresponding structure 
function. Equation (|32|) is correct in the limit of geometric optics, when the opacity variations occur on a scale that 
is large compared to D a b s /fc, which is likely to be the case in practice. However, Appendix IU1 also contains a more 
general expression for the mean square visibility when wave optics are important. 

The physical significance of Eq. (|32[1 is that the opacity variations increase the amount of fine-scale structure in the 
image of the absorbed source. This is evident through the effect of the quantity exp{[3C T (0) — D T (r)]/2}. The term 
involving C r (0) boosts the mean square visibility of the image by an amount that is independent of baseline length, 
r. However, the influence of the D T (r) partially offsets this effect, so that there is an overall enhancement in the mean 
square visibility by a factor of exp[C T (0)/2] for very large baselines. 

The structure function of opacity variations is directly related to the power spectrum of optical depth variations, 
P t (k), by a Fourier transform: 



D T {r) = 2 J d 2 K, [1 - exp(iK • r)]P T («). 



(33) 



The power spectrum of optical depth variations is well known for some spectral lines. For instance, recent studies 
of the Milky Way, Large and Small Magellanic clouds, and nearby galaxies (Braun 1999, Lazarian & Pogosyan 2000) 
show that the power spectrum of HI emission fluctuations 2 follows a power law. Green (1993) finds the index of the 
power law to be in the range n — 2 — 3, Deshpande, Dwarakanath & Goss (2000) find n ~ 2.75. Lazarian & Pogosyan 
(2000) show that the index of the power law found depends on the spectral resolution of the observations, and on 
whether any averaging has been performed over spectral channels. These authors tentatively conclude that the density 
fluctuations follow a Kolmogorov (n — 11/3) power law. Observations of HI absorption also suggest that the structure 
follows a power law, even out to scales comparable to the size of the absorbing galaxy (e.g. Braun 1999). Motivated 
by this observational evidence, we consider optical depth variations from a power spectrum of the form 



P t {k) = Q \k\ n exp 



K 



2/^ n 



(34) 



where q m ax = ^in 1S * ;rie high-wavenumber cutoff. In a turbulent medium, this cutoff physically corresponds to an 
inner dissipation scale, Z m i n , associated with the turbulent cascade. The structure function for n < 4 and r 3> ? m in 
takes the form 



D T (r) 



n— 2 



with r T = 



2 2 - n nQ - 



r (fJ 



l/(2-n) 



(35) 



2 As we refer to fluctuations in the plane of the sky, the argument of the power spectrum is a two-dimensional vector. This 
power spectrum depends on the manner in which the averaging over radial distance (or radial velocity) has been performed. 
The inversion from the full three-dimensional HI density field to the power spectrum of optical depth variations is discussed at 
length in Lazarian & Pogosyan (2000). 
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It is useful to regard the quantity r T as the length scale in the absorbing medium over which the root-mean-square 
optical depth fluctuation is unity. The optical depth structure function contains information on the distribution of 
opacity structure in the absorbing (or emitting) medium, and is a crucial ingredient in determining the image structure 
of any background source absorbed by such a medium. 



3.2.1. Scintillation characteristics of a source absorbed or amplified by a random medium 

As an application of the foregoing formalism, let us consider the practical case of variability from a spectral line in 
which the power spectrum of the inhomogeneous absorbing medium follows a power law. 

Recall that the response of image structure to interstellar scintillation is described in Sect. 2.2, from which we 
deduce that the mean square of the intensity variations is 



(AS 2 ) = 87rr e 2 A 2 SMe- 2M+3CM0 )/ 2 



V sr 



t) 



k 13 sin 2 



K 2 L 

~2k 



exp 



-Da 



kL 
~2k 



1 ( kL 
2M2* 



(36) 



In order to elucidate the properties of spectral line variability from a random absorbing medium, let us consider 
the simple case of absorption in front of a source of flux density S src and with angular diameter, # src , larger than or 
comparable to the refractive angular scale. Mathematically, this means that the source visibility declines for baselines 



at r smaller than (k0 STC ) 



< r. 



We assume that the scintillation occurs in the regime of strong scattering where 



y/L/k > raiff. Thus the term exp[— D$\ appearing in Eq. (|3 fc>|) would always cut the integrand off before the sine- 
squared term begins to vary, so the sine-squared term can be approximated by ft 4 7p/4. In practice, the upper bound 
of the integration is determined by whichever of the source or the absorbing medium cuts off the integrand at the 
lowest wavenumber. We further assume that the lower bound of the integration is much smaller than either of the two 
possible upper bounds. The resulting variation in the line flux density is 



(AS 2 



S'i 



rcf 1 



-2(t)+3CV(0)/2 



(fe) 
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(37) 



where 



S STC 2 



/3-3 



{13-2) 



IW2) / r diff 



r(2 - (3/2) V r F 



8-2/3 



rdiff < r F , 



(38) 



is the mean square fluctuation amplitude that a point source of flux density S src undergoing refractive scintillation 
would exhibit (see Fig. 3) . One clearly sees that the amplitude of the scintillation is enhanced by the additional presence 
of optical depth variations. Even when the source size is large (i.e. when 8 SIC > l/kr T ), the spectral line variations still 
differ from those of a continuum source of the same size; a continuum source would exhibit mean square flux density 
variations of amplitude S 2 of (9 src /9 SIC )^~ 6 , but the spectral variations are modified by the factor exp[— 2(t) +3O t (0)/2], 
which indicates the additional contribution of optical depth variations. The case r < (fc# src ) -1 is relevant when the 
optical depth variations occur on exceptionally short scales. For instance, for a source size src = 1 mas this case is 
relevant when r T < 3 x 10 7 A. We expect this condition to be rarely realised in practice in an absorbing medium, 
although it may be relevant to some masing systems. 

The root-mean-square flux density deviation is parametrized in terms of the optical depth variance, a 2 , and 
scattering measure, SM, as follows 



(AS 2 ) 1 / 2 



= p -{r)+3<r 2 T /i 
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(39) 



r> < {k9 sl - c y 



For scattering by Kolmogorov turbulence, corresponding to (3 — 11/3, these scaling relations reduce to 
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Fig. 3. The amplitude of refractive moduluation exhibited by a point source of flux density S src for /? = 11/3 and for 
a number of scattering measures and distances to the scattering medium, as given by Eq. lj3"%|l . 

4. Application to the variable absorption towards PKS 1127— 145 

We now apply the foregoing theory to observations of a variable absorption system. PKS 1127— 145 is a 5 — 7Jy, 
z = 1.187 quasar which exhibits continuum and HPspectral variability on a time scale of t ~ 2 days (Kanekar & 
Chengalur 2001). The intervening system absorbs 0.5 Jy, and is located at a redshift z a b s = 0.3127. After shifting all 
HI spectra to the same continuum flux level, differences of f=a 35 mJy were found in the strengths of the absorption 
lines from one observation to the next, showing that the strength of the lines varied relative to the continuum. The 
fluctuations are reported to be consistent with variations only in the depths of the components, and not with changes 
in their positions or widths. 

Kanekar & Chengalur consider HI variability due to motion of the line of sight towards a hypothetical superluminally 
ejected subcomponent of the background source. They argued that the large, rapid continuum variability on inter-day 
time scales requires S s» (3 app ~ 100. The time scale of changes in the HI spectrum implies a length scale over which 
the HI varies significantly in the intervening system. Assuming that the temporal variations are due to changes in the 
line of sight to a superluminal background source, variability on a time scale of t implies substantial changes in the 
HI structure on the angular scale of (see i|A.5|) in Appendix A) (for fi = 1 and H = 65km s -1 Mpc -1 ) 



At the redshift of the absorber, 1.7/xas corresponds to a transverse scale of 7.7 x 10 _4 pc. For the ejection speed 
argued by Kanekar & Chengalur and the observed variability time scale, this implies HI structure on an angular 
scale of 0.34 mas. This angular scale, if attributed to background source motion, is sufficiently small that interstellar 
scintillation must have at least a significant influence on the HI spectral variability. 

One can also argue that it is not feasible to explain the continuum variability as intrinsic to the source. For 
observations at the frequency of the HI absorption (i> bs — 1081 MHz), the variability, if intrinsic to the source, implies 
a brightness temperature 



The brightness temperature implied by S„(f bs) variations on a two day time scale requires a Doppler boosting factor 
in excess of S ~ 700 for the brightness temperature not to violate the inverse Compton limit in the rest frame of the 
emitting plasma. This boosting factor is unacceptably high (Begelman, Rees & Sikora 1994), and it is not feasible to 
explain such variability as intrinsic to the source. Moreover, such variations imply an observed source size 



Even if the Doppler boosting factor were an order of magnitude less than that implied by the brightness temperature, 
the variations require a sub micro-arcsecond angular size. The source would necessarily undergo scintillation, and 
would even be expected to exhibit 100% modulated intensity variations due to diffractive scintillation, which is not 
observed. We therefore attribute the continuum variability entirely to interstellar scintillation and consider a model 




(41) 




(42) 




(43) 
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based on scintillation variability alone. The values derived above are by no means unique. It should be remembered 
that the specific combination of source size, absorber size and the degree of absorption associated with the compact 
cloud are but one of a large number of possible combinations. 

In the scintillation model the continuum variations imply a source size 

"•■- 35 fe)te)(a^)""- < 44 » 

for observations at the absorption wavelength, A 0.2773 m. We apply a single cloud absorption model (see Sect. 3.1.2) 
to observations of this system to explain the amplitude of the variations observed. Let us suppose # src = 0.20 mas 
(e.g. i var = 2 days, v = 30km/s and L — 350 pc), and take the size of the absorber to be comparable to this 
at # a bs = O.lOmas, giving 9 cross — 0.155 mas and c s — 0.089mas. The flux density associated with the compact 
component in the 5 — 7 Jy source is unknown; we take the flux density of the compact component in the source to 
be Ssrc = 2.5 Jy. The absorbed flux density is observed to be w 0.5 Jy, but if only 2.5 Jy of the 5Jy source is due 
to absorption against a compact (scintillating) component one has Sabs = 0.25 Jy. For a single absorber, the given 
angular sizes and absorbed flux density imply A — 5 f a b s /S'src(^src/^abs) 2 = 0.4. The variance in spectral line intensity 
fluctuations that are independent of the continuum fluctuations is (using Eq. 1|11|0 

(A W^) 2 ) 1 / 2 = 0.070 [ W ^ m ) ^ (^-) ' Jy- (45) 

For scattering parameters typical of lines of sight off the Galactic plane, Cf^ — 10~ 4 m~ 20 / 3 and L = 0.5 kpc, the 
continuum level is expected to exhibit root-mean-square intensity variations of 450 mJy and spectral line variations of 
170 mjy The HI line, once the continuum fluctuations are subtracted, is predicted to exhibit independent root mean 
square fluctuations of 70 mJy. Thus even the simplest scintillation model is easily capable of explaining the observed 
35mJy variations in the HI line relative to the continuum in PKS 1127—145. 



5. Discussion 

Scintillation can induce variability in absorption lines relative to the continuum flux density and to other absorption 
lines. This is a natural consequence of the fact that images of the source on and off the absorption frequency differ, 
and that scintillation is sensitive to these differences in image structure. 

Scintillation effects are important whenever the source 'image' contains significant structure on angular scales 
comparable to or below the scale of refractive scintillation. Variations are therefore possible even when the back- 
ground source is not sufficiently compact to exhibit continuum variability: one only requires absorbing structure to 
be distributed on a sufficiently small scale. The angular scale probed by scintillation is typically ~ 0.1 — lmas for 
observations of HI in the local universe, which corresponds to scales from a few tens of milliparsecs to a few parsecs for 
the range of distances to extragalactic absorbers. Studies of HI absorption in our own Galaxy show that HI variations 
certainly exist on these scales (e.g. Dieter et al. 1976, Diamond et al. 1989), and there is every reason to expect similar 
variations in other absorbing galaxies. However, both the size of the background source and the manner in which 
absorbing material is distributed across it influence the amplitude of any potential spectral variability. 

Spectral line variability is expected to be more prevalent in systems absorbed against compact rather than extended 
sources. This is because the intensity variations due to structure imposed by clouds distributed randomly against a 
large background source add together independently, and reduce the amplitude of the modulations relative to the 
total absorbed flux density. To see this, consider the variability due to a collection of N clouds absorbing radiation 
against an extended background source. The power spectrum of intensity scintillations across the observer's plane 
is proportional to the power spectrum of the brightness distribution. The amplitude of the brightness distribution 
scales only linearly (rather than quadratically) with the number of absorbers if they are distributed randomly and 
sparsely across the background source. Thus, for a fixed absorbing flux density the root mean square amplitude of the 
modulations scales as TV" 1 / 2 . This is of particular relevance to very large background objects, where the number of 
foreground absorbers can be correspondingly large. 

When the background source is compact its structure acts like a high-pass filter on the structure of the absorbing 
material, rendering spectral variability more likely. In the case of absorption from discrete clouds, this is because the 
number of independent clouds contributing to the absorption is reduced. A background source which is itself compact 
enough to exhibit interstellar scintillation is guaranteed to exhibit spectral line variability relative to the continuum 
variations provided there is some variation in the optical depth across the source. 

Scintillation is only capable of modulating the intensities of existing features; it does not alter the frequencies of 
the spectral lines. However, this does not preclude variability in the shapes of line profiles. This is expected if the 
angular distribution of absorbing material across the background source varies across the line profile, and is particularly 
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relevant to observations in which contributions from individual absorbing features overlap spectrally. This would be 
the case if many distinct clouds contribute to the absorption at a single frequency. In essence, profile variability can 
occur because the apparent brightness distribution varies within the line profile and causes certain parts of the line 
profile to be enhanced relative to others; the scattering medium sees different source 'images' at different frequencies. 
Thus, while scintillation can not shift the absorption frequency, it can alter the line profile by enhancing contributions 
from one frequency relative to another. 

When a spectral feature contains contributions from many absorbing clouds, the outer edges of the line profile are 
expected to exhibit greater fractional variations relative to the line centre. The edge of the line profile is more likely to 
contain very compact structure, due to small amounts of material at extrema of the velocity distribution. Indeed, such 
behaviour may have already been observed at the edge of the Hi-absorption line profile in PKS 1127—145 by Kanekar 
& Chengalur (2001). 

Different absorption lines from within the same system arc expected to show differing degrees of variability, ac- 
cording to the amount of fine scale structure associated with each absorption line. Parts of the ISM at different 
temperatures in an intervening absorbing system arc likely to possess different angular distributions. A colder, more 
densely clumped medium with more power on smaller angular scales would exhibit more variability. Such a naive 
model would predict narrow lines to exhibit more variability relative to the broader lines from slightly warmer parts 
of the medium. 

The response of scintillation to a power law spectrum of absorption (emission) fluctuations is determined by the 
variance and scale length of the optical depth fluctuations in the medium. For a background absorbing source larger 
than or comparable to the refractive scale, the presence of the optical depth fluctuations enhances the mean square 
amplitude fluctuations by a factor exp(3(r 2 )/2). This enhancement is due to the extra structure imprinted on the 
background source by the absorbing medium. 

Scintillation- induced spectral variability is more important for absorbers at high redshifts. Absorbers of a fixed 
spatial scale appear smaller at high redshifts. The apparent angular diameter of structure of a fixed linear scale, So, 
decreases with redshift as s /Da(z) for angular diameter distance Da(z). As the angular diameter distance increases 
until at least z ^ 1 (depending on the cosmology assumed), so does the apparent source size decrease. A second 
effect also enhances spectral variability at high redshift. The angular scale probed by the refractive scintillations 
increases with wavelength as A 2 2 . Scintillation effects are most important when this scale exceeds the angular scale 
of the absorption structure. It is easier for absorption which occurs at higher redshifts, and thus at longer apparent 
wavelengths, to satisfy this criterion. 

Several other possible causes of spectral line variability have been suggested. Lewis & Ibata (2003) discuss the 
contribution of gravitational lensing to absorption line variability. Like scintillation, variability due to gravitational 
lensing requires the presence of small-scale structure in the absorbing material. However, the presence of lensing 
material is problematic: it is not clear whether the lensing optical depth in intervening systems is sufficiently large to 
reproduce the continuous sort of variability that is reported. On the other hand, if lensing effects are important then 
it is likely that scintillation variability will also contribute strongly to the variability, as both mechanisms depend on 
the presence of small scale structure in the absorbing medium. 

Spatial variations in the absorbing structure coupled with motions in the background source can also cause line 
variability. Both Kanekar & Chengalur (2001) and the present paper discuss the possibility of source motions giving rise 
to HI variability in PKS 1127—145, while Briggs (1983) considers this in relation to the HI variations observed in AO 
235+164. For the model considered by Briggs (1983), the optical depth variations are not independent between spectral 
lines. Measurements of any two lines are sufficient to predict the variations of the other lines observed. Application 
of this model to AO 235+164 shows that background source motions can reproduce the spectral variability observed. 
However, we note that the presence of correlations between spectral line variations does not preclude scintillation as 
a viable explanation, as scintillation causes correlated variability between sets of spectral lines whenever there is a 
correlation in the spatial distribution of the absorbing material. Indeed, if the structure were all located within a small 
enough angle (i.e. within 6 Te f), the variability would be completely correlated. 

6. Conclusions 

We have presented a theory for the effect of interstellar scintillation on the variability of HI absorption lines from 
extended (i.e. non point-like) extragalactic sources. The model predicts modulations of the line strengths relative to 
the continuum and relative to other lines. It does not permit shifts in the frequency of the lines. However, apparent 
shifts in the line frequency can occur if the angular distribution of the absorbing structure differs with frequency within 
the line. 

Scintillation is most sensitive to structure on scales less than the refractive scale of the scintillation, typically ~ 
1 mas. Line variability arises when intervening absorption imposes structure in the source image that varies significantly 
with frequency on this angular scale. 
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The theory presented here applies equally to any absorption or emission line in the radio band provided that its 
distribution is sufficiently compact for scintillation to be important. In this respect, the foregoing theory is also useful 
in the interpretation of fast maser variability (Greenhill et al. 1997). 

Other possible explanations of spectral line variability include gravitational lensing effects and motions in the 
position of the background source relative to the absorbing material. However, the presence of small scale structure in 
the absorbing material in both cases means that scintillation effects are likely to play an important, if not dominant 
role in the variability. Scintillation is a viable explanation for the cases of spectral variability hitherto observed at 
radio wavelengths, particularly for observations of the HI line. 

As an application of the theory, we have examined the HI variability reported along the line of sight to PKS 1127—47, 
in which spectral lines vary relative to one another and to the continuum. A simple model accounts for all the observed 
inter-day spectral variability. At this frequency, scintillation can account for variations on time scales of days to weeks. 
We also argue that scintillation must be important in this source even as an explanation of the continuum variability. 

Acknowledgements. I thank Nissim Kanekar for many useful discussions and for reading the manuscript critically and Avinash 
Desphande for reviewing the manuscript and pointing out a conceptual error in the original manuscript. 



Appendix A: Time scales for superluminal motion of a background source relative to an intervening 
absorbing medium 

Often the emergence of a new compact background component is associated with the relativistic ejection of material 
from a background quasar. The apparent transverse velocity of the component often exceeds c, and one expects the 
point at which the line of sight intersects an intervening absorption system to change with time. Here we compute the 
time scale of variability due to such motion. 

Suppose an observer at Earth measures a change in angle /i in the direction to the background source over a time 
interval. The change in angle corresponds to a distance x src = (j,D A (z SIC ) at the source, where Da(z sic ) is the angular 
diameter distance, given by 



2c f n z + (n-2)(VT+m-i) \ . 



H \ ft 2 (l + z) 2 

The rate of change of the angle fi (i.e. the proper motion) measured by an observer on Earth can be related to the 
velocity 

dt E dt E Da{z s ?c) 



dtE dt STC \Da(z stc ) 

Vsrc 

(1 + Z sxc ) D A (z SIC )' 



(A.2) 



We have used the fact that frequencies in the source frame and Earth frame are related by fE/^src = dt SIC /dtE = 
(1 + z)^ 1 -. an event of duration At E measured on Earth corresponds to a duration A*e/(1 + z) in the frame of the 
emitting system. A source expanding at v S rc = (3c as measured at the redshift of the source has a proper motion 

A = 777TT y (A-3) 

D A {1 + z src ) 

The proper motion is determined by the apparent speed of the material, which may exceed c by a factor of tens, and 
by the redshift. We can also use Eq. (|A.2|I to deduce the proper velocity, u a bs, at which the line of sight moves across 
an intervening HI absorbing screen, located at redshift z a b s 

Wabs = (lD A {z a bs)(l + Zabs)- (A. 4) 

The speed at which the line of sight moves across the HI absorbing screen is thus 

Vabs = /3app C — — r — . (A.5) 

L>a{Zsic) \ 1 + Z SIC J 
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Appendix B: Average power spectrum for a large number of absorbers in front of a compact source 

If there is a large number of absorbers covering the background source it is appropriate to average over their positions. 
The average over AOi is performed by assuming they are randomly distributed in space out to some angular distance, 
#max with a probability density function 



p(A6) 



(2#max) 2 , \A9 X \ < 

■/max i 

\A6y\ < 

0, otherwise 



The average over positions is performed by integrating over the probability distribution 
" ,ax dA6 x dA6 v 



46> 2 

max 



exp [i(a x A9 x 



1666*2 



■ exp 



a y Ae y )-b(A0 2 x + A%)] 
2b0 m 



-lb 



erf - 



2y/b 



erf 



2b6 n 



2Vb 



where 
(fix, a y ) 
b 



9 2 

y src 



9 2 

src 



erf 



(k x , K y) 



2b9„ 



ia,, 



2Vb 



erf 



266L 



2Vb 



or 



(B.l) 



(B.2) 

(B.3) 
(B.4) 



Appendix C: The average mean square visibility due to a random absorbing or emitting medium 



The wavefield incident upon a point r on a plane at a distance D from a thin screen of absorbing material with phase 
variations cf> can be written in terms of the Fresnel-Kirchoff integral (see, e.g., Born & Wolf 1965, Goodman & Narayan 
1989) 



u{r) 



JkD 



2nir F 



J <i 2 xu src (x) exp 



.(x-r) 2 
'^4~ 



(C.l) 



where u src (x') is the wavefield of the background source incident upon the absorbing medium and the Fresnel scale 
is given by rp = y^D/k. The background source is assumed to be located at infinity, so that its wavefronts incident 
upon the absorbing plane are planar. 

Here (j) is regarded as a complex quantity. Since we are interested only in opacity variations, we henceforth set 
Re((/>) = and consider only the imaginary contribution to the phase, lm(0) = rj = t/2. With this specification, Eq. 
IjC.ip is a solution to the wave equation for a thin screen of absorbing material in the absence of refraction. 

Equation (|C.1|1 incorporates wave effects that are normally neglected when considering the effect of optical depth 
variations. For instance, if one points a telescope toward a particular point in an absorbing medium, one expects to 
measure an optical depth associated with that particular point in the medium. However, although the optical depth 
t is a well-defined quantity, Eq. I|C.1|) implies that the wavefield one measures contains the contributions r averaged 
over some region in the absorbing material, with the weighting function given by the first exponential term appearing 
in the integrand. This effect is important when the optical depth variations occur on a scale that is small compared to 
the length scale r-p = ^D/k. Equation l|C.l|l reduces to the familiar definition of optical depth u(x) = u src (x)e~ T ( x )/ 2 
in the limit of geometric optics, which corresponds to A — > 0. However, we choose to work within the framework of 
physical optics in the present case because its generality allows one to consider the effect of absorbing structures that 
may be small compared to rp. This case may be relevant for absorption by some extragalactic systems where D is 
large and the medium may be sufficiently inhomogeneous. 

The mean square visibility of radiation from a distant absorbed background source which has propagated through 
the random absorbing medium is 



(V ahs (r')V: hs (r')) = (u{v + v'/2)u*{v-v'/2)u{v-v'/2)u*(v + v'/2)) 

= 7^ — 2T7 / rf 2 Xid 2 X2d 2 x 3 d 2 x 4 (u src (x 1 )u* rc (x2)u src (x 3 )u* rc (x4)) (e 
(27rrf-) 4 J \ 

. (xi - r - r'/2) 2 - (x 2 - r + r'/2) 2 + (x 3 - r + r'/2) 2 - (x 4 - r - r'/2) 2 



5 -J)(xi)-J)(x 2 )-?)(x3)-?;(x4) 



x exp 



2r\ 



(C.2) 
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e -'/(xi)-i)( x 2)-'/(x3)-r/(x4) \ 



I'm' a spatially incoherent background source, the average over the source wavefield appearing in the integral reduces 
to a product of visibilities, V(xi — x 2 )F*(x 3 — x 4 ). We use the result (exp[— x}) = cxp[— (x) + (x 2 }/2] for a normally 
distributed variate x to reduce the average over opacity variations: 

exp { - 4(77) + - [120,(0) - A,(xi - x 2 ) - A,(xi - x 3 ) - D v {xi - x 4 ) 

-A,(x 2 - x 3 ) - A,(x 2 - x 4 ) - A,(x 3 - x 4 )] |, (C.3) 

where C n {r) = ([rj(r + r') — (77)] [rj(r') — (77)]) is the autocovariance of opacity variations and D v (r) — 2[C V (0) — C v (r)] 
is its corresponding structure function. 
Making the change of variables 

R- = j [xi + x 2 + x 3 + x 4 ] 

! r 1 
ri = - [xi + x 2 - x 3 - x 4 J 

! r 1 
r 2 = - [xi - x 2 - x 3 + x 4 J 

p = xi - x 2 +x 3 - x 4 , (C.4) 

the mean square visibility reduces to 

(u(r + r'/2)u*(r - r'/2)u(r + r'/2)u*(r - r'/2)) 
1 



d Rd pd Yi d r 2 exp 



>i • (r 2 - r') + p • R - p ■ r 
1 ^ 



V Br c(r 2 +p/2)V r Brc (/o/2-r 2 ) 



" (2^-2)4 

exp { - 4(77) + 1 [120,(0) - A,(r 2 + p/2) - D v {r x + r 2 ) - D,(ri + p/2) - D^(r x - p/2) 

-D v {vx - r 2 ) - D n (p/2 - r 2 )] }. (C.5) 
The integral over R yields a factor 5(p)(27rr|) 2 , and the subsequent integration over p eliminates r, yielding, 



(u(r + r'/2)u*(r - r'/2)w(r + r'/2)u*(r - r'/2)) 

. (r 2 - r') ■ ri 



(27rr|) 2 



d ri d r 2 exp 



^src(r 2 )K; c (r 2 ) 



exp { - 4(7;) + i [120, (0) - 2D v (n) - 2A,(r 2 ) - ^(ri + r 2 ) - D„(n - r 2 )] }. (C.6) 

In the geometric optics limit wave effects are removed by taking the limit A — > 0, which corresponds to the limit 
tf — > 0. In this case one can show, either from Eq. IC.ljl or directly from Eq. (|C.6jl . that the mean square visibility is 

(|Kb s (r')| 2 ) = |^src(r)| 2 exp [-4(77) + 6C„(0) - 2D v (r')} . (C.7) 
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